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O \ Abstract 

The scattering process on multiloop infinite p + 1-valent graphs (generalized 
^ , trees) is studied. These graphs are discrete spaces being quotients of the uniform 

^ I tree over free acting discrete subgroups of the projective group PGL{2, Qp). As the 

^ \ homogeneous spaces, they are, in fact, identical to p-adic multiloop surfaces. The 

Ihara-Selberg L-function is associated with the finite subgraph — the reduced graph 
Q>^ ■ containing all loops of the generalized tree. We study the spectral problem on these 

' graphs, for which we introduce the notion of spherical functions — eigenfunctions of a 

. discrete Laplace operator acting on the graph. We define the S'-matrix and prove its 

unitarity. We present a proof of the Hashimoto-Bass theorem expressing L-function 
, of any finite (reduced) graph via determinant of a local operator A(ii) acting on 

Ch \ this graph and relate the S'-matrix determinant to this L-function thus obtaining 

^ I the analogue of the Selberg trace formula. The discrete spectrum points are also 

determined and classified by the L-function. Numerous examples of L-function 
calculations are presented. 
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Introduction 



Two fundamental results relating the L-functions (zeta functions) and Laplacian determi- 
nants have been obtained respectively by Selberg |p2| (developed by Lax and Phillips , 



L. D. Faddeev and B. S. Pavlov |]T2[, A. B. Venkov |3^, and others) and Ihara pO| , 
Hashimoto \TE\, and H. Bass 0. The first result — the celebrated Selberg trace formula — 
concerns the zeta functions and Laplacian determinants of compact (or noncompact with 
finite area) in general multidimensional (D-dimensional) manifolds of constant negative 
curvature, which can be obtained from the D-dimensional upper half-space with the hy- 
perbolic metric by factoring it over a discrete acting subgroup of the symmetry group of 
this space. We briefly expose these results in order to compare them with the second class 
of results concerning the Ihara-Selberg L-functions on finite graphs. 

The zeta function associated with a finite graph was introduced by Ihara who 
proved the structure theorem for torsion-free discrete cocompact subgroups P of the group 
PGL{2, Kp); Kp being a p-adic number field or a field of power series over a finite field. 
Then, P is a free group and the associated zeta function can be constructed as follows. An 
element 7 G P, 7 7^ 1, is a primitive element if it generates its centralizer in P. Define the 
length l{'~f) = ||Ai/A2||p of the element where || ■ ||p is the normalized valuation in Kp and 
Ai, A2 G Kp are the eigenvalues of 7. Let V{T) be the set of conjugate classes of primitive 
elements of P. Then, the Ihara zeta function is 



-revir) 



(0.1) 



This definition was extended to finite-dimensional representations of P in [^]; then, in- 
stead of (10.11), we have 



Zriu,p)= n det(/i-«'Wp(7))"\ 

7GP(r) 



(0.2) 



where ^(7) is the character of the given representation and Ii is the unit operator. (Orig- 
inal formula (|0| ) corresponds to p{-y) = 1.) 

Although zeta function (|0.2| ) is an infinite product, it was proved to be a rational 
function. The Ihara theorem states that 



Zy{u, p) = (1 - u^y det(/o - Au + pu^ 



(0.3) 



where s is an integer and A is a matrix acting in a finite-dimensional linear space. Us- 
ing the rationality formula one can counts the number of conjugate classes of primitive 
elements of P. 

1. In 1^ scattering processes on infinite graphs representing spaces of a constant 
negative curvature were studied. Such are uniform, or univalent graphs, for which a 
natural number p > 1 exists such that all vertices of the graph are incident exactly to 
p + 1 ends of edges. For prime p, such graphs can be interpreted as homogeneous spaces 
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for the p-adic projective group PGL{2,Qp) factorized, first, by its maximal compact 
subgroup PGL{2, Tip) and, second, by some discrete, free-acting (Schottky) group V. 
The homogeneous space D ~ PGL{2,Qp)/ PGL{2,Zp) is a uniform infinite tree graph. 
If we now factorize the tree D by a discrete freely acting finitely generated subgroup 
Tg C PGL{2, Qp), where g is the number of the generating elements, then the obtained 
graph Tg = D/Tg is again a univalent graph containing g loops and having the tree D as 
the universal covering. 

We present a proper analogue of the Selberg trace formula |^ for such discrete sur- 



faces. These graphs sometimes are referred to as multiloop p-adic surfaces because it has 
been shown in [^, that for p the prime, the discrete Laplacian action on these graphs 
yields the proper scattering amplitudes of the p-adic string. 

Note, however, that most of the calculations were performed for trivial (Abelian) 
representation of the group F in the p-adic string approach. 

A natural distance \Px,y\ over a path P^^y connecting the points of a graph is merely 
the number of edges entering (with multiplicities) this path. One can consider linear 
spaces of functions Cq and Ci depending, respectively, on vertices Xi and oriented edges 
Cj of the graph T.. A Laplace operator A acts on the space Cq in a standard way, 

A/(^)= E f{x.)-{p + l)f{x). 

I Px , X , j I — 1 

It is useful to segregate from the graph T its finite "closed" part — the connected reduced 
graph F containing all internal loops. The valences of F vertices can be arbitrary (< p+ 1). 
In 1^ , the string theory for such graphs was developed and the proper p-adic counterparts 
of all crucial ingredients of the open string theory, such as prime forms, Schottky groups, 
etc. were found for the corresponding scattering amplitudes. 

2. For the ordinary closed (or, more general, finite area) Riemann surfaces of con- 
stant (negative) curvature the Selberg trace formula, which establishes an explicit relation 
between determinants of the Laplace operators and zeta functions (or, Ihara-Selberg L- 
functions), is indeed close in form to ( p.3|) . The distance ^(7) is now the length of a closed 
geodesies in the constant negative curvature metric. However, no proper analogue of the 
Selberg formula is known for the essentially noncompact case where the spectrum of the 
Laplacian contains a continuous part responsible for the scattering processes. Note, how- 



ever, that zeta function (|0?T|) is still well defined for such noncompact surfaces; it is the 
Laplacian spectrum definition that hinders the progress in this direction. One may hope 
that the graph description of moduli spaces of complex curves with holes p3, ITl, B] related 



to the geodesic description may help in finding analogous statements for the rationality 
in the Riemann surface case. We discuss this approach in Sec. 6. 

3. We adopt the definition that zeta functions take values in the field kp (on the 
definition of the corresponding noncommutative determinants, see 0), while L-functions 
are assumed to take values in the complex number field C. 



In a series of elegant papers by K.-I. Hashimoto [|T^] and H. Bass on zeta and 
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L-functions, the Ihara formula was generalized to the case of an arbitrary (not necessarily 
univalent) finite graph. In the most general case, it has the form 



L{u,p) = det-\li-uTi), (0.4) 

where the new operator Ti is the operator acting on the space of functions on oriented 
edges (simplexes of dimension one), which corresponds to a translation along edges of the 
graph. 

Relation ( p.4|) is suitable for both zeta and L-functions; the difference is only in the 
definition of determinants (commutative or noncommutative) involved. In this paper, 
we consider only the case of L-functions, which seems to be physically meaningful as 
pertaining to the spectral and scattering problems on graphs with potentials. 

Formula (p.4|) is amazingly universal; it expresses the rationality condition, which holds 
for the zeta functions in a host of cases related to discrete dynamics — one may generalize 



it to the case of Hermitian (non-Abelian), Schrodinger type potentials on edges |26], to 
the case of nonzero torsion (potentials depending on "angles" along the path, i.e., on the 
pairs of consecutive edges (64,64+1)) [||, etc. However, in most of these cases, we cannot 
find a bridge from ( |(J.4| ) to the Laplacian determinants on nonuniform graphs. Only when 
the potential Uft on (oriented) edges satisfies the unitarity condition, i.e., f/^ = U^^ for 
all pairs of oppositely oriented edges, we can express the L-function in terms of Laplacian 
type determinants on arbitrary (not necessarily univalent) graphs, 

L{u,p) = {1 - uy det-\l - uMi + u'^Q) = {1 - det-^A{u), (0.5) 

where Mi is the above operator of summation over all neighbors (with probable potential 
dependence) and the new in comparison with (p.2| ) element is the operator Q, which 
counts the number of these neighbors: Qx = qx if the vertex x eT has q + 1 neighbors. 
In the trivial representation case, s = \V\ — \E\, where \V\ and \E\ are total number of 
vertices and (nonoriented) edges of the graph (in a non-Abelian case, s is the differences of 

2 1 VI \E\ 

total (half) dimensions of linear spaces vj^ and , where elements of the representation 
v_A of the group or algebra A dwell on (now oriented) edges and vertices of the graph. 
However, the operator A{u) becomes the Laplacian only for u = 1; these operators even 
do not commute at different u. (For the detailed description of these results and their 
different applications, see [0.) 

4. The graphs T (p-adic multiloop surfaces) are noncompact in general, which differs 
them both from the closed Riemann surfaces and with finite graphs. However, it is possible 
to use the spherical function technique and describe the problem in physical terms of the 
scattering theory on graphs in this case. In the original setting, spherical functions 
are eigenf unctions of the Laplace-Beltrami operator that depend only on the distance to 
a given point xq (the center). Because the Laplacian is the second-order operator, we 
always have two branches of the solution (at a distant point) proportional to a'^ and al 
where d is the distance to the center. Resolving the eigenvalue problem at the central 
point we fix the ratio of coefficients a+ and a_ standing by these two branches. If we 
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choose "^(^^o) = 1, then a+ and a_ become Harish- Chandra coefficients and c = a+/a_ is a 
scattering amphtude of the s-wave. Spherical functions have been found for the scattering 
on a quantum hyperplane |]35| as well as for the scattering on p-adic hyperbolic plane |T5 . 



The S'-matrices obtained are closely related to the partition functions of the XX Z model 
and a number of nice but still mysterious relations between them have been obtained by 
Freund and Zabrodin 111 



Investigating the spectral theory on graphs has already rich history. Results obtained 
were related to studying eigenvalue problem on finite uniform {i.e., univalent) graphs. 
Here, it was proved that a deep relationship exists between the modular forms on Te- 
ichmiiller spaces and finite graphs, namely, the so-called Ramanujan graphs. We do not 
discuss this interesting approach in this paper and instead refer the reader to the mono- 
graph ||31| and papers [^, |3^ . 



In [^, an analogue of a spherical function for the multiloop graph was introduced. 
The main idea is the following. Considering the spectral problem Lip = Xip, L being 
the Laplacian, we note that any linear superposition of spherical functions with the same 
eigenvalue A but different scattering centers is again an eigenfunction. An eigenfunction of 
the Laplace operator on the factorized tree T = X/Tg corresponds to a source distribution 
function, s{x), on the tree X such that for every 7 G and x E X,we have s^jx) = s{x). 
Then the whole eigenfunction is periodic under the action of F^. Moreover, we choose a 
finite domain (a reduced graph) F G T and consider only s{x) such that supp s{x) C F. 
Inside T there is a unique minimal finite connected subgraph containing all loops — the 
union -D(F) of invariant axes of all elements of F^ factorized over the action of the group 
Tg. This graph, D(F)/F, contains all information about the "geometrical structure" of T. 
We always assume D(T)/T C F. 

Each eigenfunction iIj{x) may be presented as a sum of retarded and advanced wave 
functions: 

^l;{x) = A,,,(«)af - A,,,(n)a?^'"(^» ^ ^+(x) - ^_(x), (0.6) 

where a+ and a_ are two fixed complex numbers depending only on the eigenvalue t of 
the Laplacian L = A(l), which acts on the whole graph T, Lip = {t — p — l)ip and on the 
initial prime number p, a± = ^ ± \J ^ ~ ot+(^- = and u{x) is the closest to x 
point of the reduced graph F (for x G F, u{x) = x). As in the central symmetric case, 
ip{x) depends only on distance to F on branches outside F. 

Now we can define a scattering matrix S for such system. As a basis we choose 
functions Al^^^{u) and Al.^^{u) that are nonzero only for some (ith) point of F. Then, 
determining the asymptotic vectors ipf that behave as a^. when going along ith branch 
and are zero otherwise, we can asymptotically expand ip in the sum of basis vectors 
ip^^ — Siajb'^j' ^ where now ia,jb are multiindices indicating components a of the 
representation of v_a on the ith branch. Then, the matrix Si^j^ has the natural sense of 
the S'-matrix. 

The determinant of the matrix S depends on the spectral parameter, contains the 
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information about spectrum, and, moreover, is directly connected with the L-function of 
the graph T. (As for eigenfunctions themselves, the Lax-Phillips approach was developed 
for their description pO| .) 

In order to find the determinant of S we impose a restriction: Aadv{u)/Aret{u) = const 
for all points u & F. Imposing this condition at all — both external and internal points of 
F, we fix an arbitrariness in the splitting of '4){x) into the advanced and retarded waves. 

Because the central object — the reduced graph F — is finite, a finite set of possible 
eigenvalues of the 5*- matrix, i.e., such constants q (the letter "c" originates from Harish- 
Chandra c- function |19|) exists. Their product C is therefore the determinant of the 
S'-matrix. We also called it the total C-function. We present the proof of the basic 
formula establishing the relationship between C and determinants of a local operator 
A(m) acting only on the reduced graph F, 

^/ X det A(a_) ,^ „, 

det5(t)= — ^ — f, (0.7) 



detA(a^^ 

where is the total dimension of linear space of functions at vertices of the reduced 
graph F, which take values in vji,. 

On the contrary, the operator A(a::p) in ( p.7|) is taken from ( |0l5|) , A(-u) = l+Qv?—uMi 
and it is determined completely in terms of the reduced graph itself, not of the whole 
graph; the only remaining dependence of the "big" graph T is contained in arguments 
and a_ of the function A(a±). Comparing ( p.7|) and (|075| ) we obtain the relation between 
Harish- Chandra total C-function and the L-function of the p-adic curve: 

^ot- ) \l — a\ j L{a- 

where tq — ri is the difference of total dimensions of spaces of functions determined on 
vertices and (nonoriented) edges; it is equal ((7 — 1) x where g is the number of loops 
(the genus) of the graph and \vy^\ the representation dimension. We also prove using the 
geometrical setting of that the S'-matrix is unitary in the scattering zone. 

From another point of view, the idea to consider a proper product of scattering co- 
efficients calculated at different p is originated from the adelic ideology. Such products 
in scattering processes was first proposed in |T^, where the product of C-functions for 
scattering on p-adic hyperplanes taken over all primes p appeared to be connected with 
the C-function of scattering on genus one modular figure considered by L. D. Faddeev 



C=[^] (^1 (0. 



and B. S. Pavlov ||T2|. The very general formulas concerning the scattering on symmet- 
rical spaces and gamma-function technique can be found in [|I7|. One could hope to find 
a proper adelic products of L-functions appearing in our approach in order to compare 
them with the ones for Riemann surfaces. 

We also discuss the eigenvalue problem as regarding to the discrete part of the Lapla- 
cian spectrum and establish two important relations concerning the S-matrix. The dis- 
crete part of spectrum may contain apart of "customary" discrete eigenvalues correspond- 
ing to exponentially decreasing in branches eigenfunctions also the so-called exceptional 
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eigenvalues corresponding to eigenfunctions that vanish identically on all branches. We 
first prove that poles of L-function correspond to normal discrete eigenvalues iff these 
poles do not cancel each other in (p.8|), i.e., if and a_ are not simultaneously the 
poles of the L-function. On the contrary, as soon as such situation takes place, i.e., there 
are such poles a+ and a_ of the L-function that = 1/p, the exceptional discrete 

spectrum appears. 

The paper is organized as follows: Section 1 contains definitions, the interpretation of 
p-adic multiloop curves as graphs and the action of the Shottky group on the initial tree 
graph. In Sec. 2, we describe the automorphic functions and potentials on these graphs 
and introduce operators acting on the generalized tree simplicial complexes. In Sec. 3, 
we describe L-functions associated with such groups, or, equivalently, with the reduced 
graphs. The theorem by Hashimoto and Bass is formulated and a proof is presented. In 
Sec. 4, we consider the spectral problem, introduce the spherical functions on multiloop 
graphs, define the corresponding S'-matrix and show that its determinant can be expressed 
as a ratio of two L-functions. We prove the unitarity of the S'-matrix and describe how to 
find discrete spectrum eigenvalues using the L-function technique. Examples for genus 1 
and 2 are presented in Sec. 5 together with a simple algorithm for calculating L-functions 
in lower genera based on the Hashimoto-Bass theorem. Eventually, in Sec. 6, we present 
a construction of graphs for describing the Teichmiiller spaces of complex curves in the 
Poincare uniformization picture and set the problem of finding a proper analogue of the 
Selberg trace formula for open Riemann surfaces. 

1 Definitions 

1.1 Graphs and trees 

Let p be a natural number and D a uniform tree graph of order p + 1, V{D) and L{D) 
the sets of its vertices and (oriented) edges. In what follows, we always consider oriented 
edges, i.e., a two-dimensional subspace corresponds to each (nonoriented) edge. For each 
two points X, y of the tree, the distance d{x, y) is equal to the length of the unique way 
connecting these two points. 

Definition 1 Let T be a graph with finite number of loops and branches (tails), V the 
set of its vertices, and L the set of (oriented) edges. The graph T is a coset of its universal 
covering tree, over the action of freely acting subgroup F of the group of motion of the 
tree, T = D/V. We denote by e and e two edges from L with the opposite orientations. 
The modulus | ■ | of a set is the cardinality (perhaps, infinite) of this set. 

Definition 2 Let an oriented path P^^y in a graph T (or D) be a (unique) sequence 
(finite or infinite) (ei,...,em) of consecutive {doCi = (iiCi-i, 1 < i < m, where do and 
di are the respective operators that project an oriented edge to the vertex it starts or 
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terminates) oriented edges without backtracking (i.e., 7^ ej_i for 1 < i < m) starting 
at the vertex x = doCi and terminating at the vertex y = diCm- The path length \Px,y\ 
is the number of edges entering the path. In the tree D, the length of the path P^ y is 
always the distance between the vertices x and y. 

The path P^^y is closed if x = y. The proper closed path is a closed path with ei 7^ 
(for a path of nonunit length). 

We say that an edge follows an edge e if starts at the vertex at which e terminates, 
i.e., doc^ = die. 

Remark 1 The set of proper closed paths with the marked (starting=terminating) points 
removed (by using the forgetting mapping) is in a one-to-one correspondence with the 
conjugate classes of the homotopy group vri(T). 

1.2 F-action on trees D 

We recall some facts about the construction of the action of finitely generated freely acting 
group r on the tree D. 

The tree graph D possesses a rich group of isometrics — the transformations preserving 
the distance on the graph. Considering a "rigid" graph, i.e., a graph with the fixed 
ordering of edges in each vertex, and imposing the condition of the ordering preservation 
under the action of the symmetry group, we reduce this huge group to the group of 
rotations and translations along the axes (infinite lines) of the tree. (It is possible to 
interpret the boundary of the tree D as a p-adic projective plane Pi(Qp) |2^. Then the 
full group of motions of the tree X is the projective group PGL{2,Qp).) A rotation 
element (an analogue of an elliptic element of the projective group) is a tree rotation 
about a vertex xq G -D through the angle 27m/ (p + 1), 1 < n < p. This is not however the 
case we consider in this paper. Another type of transformations is provided by elements 
that have no fixed points inside the tree (free-acting elements). 

Without losing the generality, we assume that the group F is a discrete free-acting 
finitely generating subgroup of PGL{2, Qp), that is, no fixed points inside the tree exist 
for all nonunit elements of this group. (Thus, it is an analogue of a Fuchsian group 
in ordinary hyperbolic geometry.) We consider in what follows only finitely generating 
groups F and assume that the group has g generating elements. Each element 7 G F, 
7 7^ 1, induces a translation of the tree as a whole along the invariant axis D{'y) C -D of the 
element. The invariant axis -D(7) is a unique infinite oriented path . . . e_ieoeie2 . . . that 
maps to itself under the action of 7: 7(6,) = ej+/ shifting in the positive direction w.r.t. 
its orientation. Here / = l{'~f) is the element length, ^(7) = ini{d{x,'y{x))\ x G V{D)}, 
where obviously the minimum is reached on the set x G D {'-/). Explicitly, an element 7 
defines a translation on the distance ^(7) along the line D{^). (We assume Z(7) < 00. )0 

^Note that no analogue of parabolic element exists in this geometrical setting. 
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A centralizer ^'(7) C F of the element 7 G F is a cyclic group. If, besides, 7 is a 
generator of Z{'y), then 7 is called a primitive element of F. It is useful to define a von 
Mango It function A (7) = /(tu), where 7 = tu™ and is a primitive element of F. 

Let us consider now the subtree -D(F), which is the union of all invariant axes of the 
elements of F. For an element 7 having the length /(7), the sequence . . .e_ 1606162 . . . 
has the structure . . . (6*162 . . . e/)(e*]^e*2 . . . e])(e*i' . . . e]') . . ., where e^"^ are copies of the edge 
Ci G T. Moreover, the sequence 6162 ... 6; must be the proper closed path in T (with 
possible repetitions, i.e., some of may have the same preimages in T). 

Let us consider the set of conjugate classes {7} of the group F: 

{7}: I U ^7^-'eFl. (1.1) 

For each element (3 G {7} and each vertex x G V{D), we have (3x = wyu'^x = y; hence, 
{u!~^y) = ■y{u!~^x). Therefore, for each element /? G {7}, there exists such c<j G F that the 
invariant axis of (3 is an image of the invariant axis of the generating element 7 under the 
action of uj~^: D{P) = uj~^D{^). Thus, there is a one-to-one correspondence between 
the conjugate classes {7} and a set of all proper closed paths in the graph T. 

We call primitive conjugate classes {w} such classes {7} that are generated by the 
primitive elements of F. Then, the proper closed path corresponding to {w} is the proper 
closed path in T that has no subperiods. 

All graphs we consider are obtained by factoring a tree D over the action of a free- 
acting symmetry group. 



1.3 Reduced graphs and branches. 

Definition 3 Let the reduced graph F be a finite, necessarily connected subgraph F G T 
containing all loops of the graph T. Its universal covering. Dp, is a subtree in D. 

For the generality sake, we do not demand the reduced graph to be inambiguously 
determined by the graph T. However, we always can segregate the minimum reduced 
graph, which is the intersection of all reduced graphs admitted by the given graph T. This 
minimum reduced graph exactly coincides with the union of all -D(7), 7 G F, factorized 
over the action of the group F. This union is a subtree -D(F) C D, but neither -D(F) 
nor F must be uniform graphs. 

The subgraph D{r)/r and, correspondingly, a reduced graph T are finite graphs 
containing exactly g loops, where g (the genus) is the number of generating elements 
of the group F. 

Because a subtree -D(F) does not coincide in general with the whole tree D, a quotient 
F\X would be an infinite graph in contrast to graphs in paper . 
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The graph T can be presented in the form T = F U B{T), where all g loops of the 
graph T are contained in the reduced graph F. Meanwhile, the minimum reduced graph 
contains no terminal points, i.e., such points that are incident to only one edge in F. The 
complement to F, B{T), is a (finite or empty) set of branches growing from the vertices 
of -F in a way to make the total valence of a vertex of T to be p + 1. We always assume 
that the number of branches growing from a vertex of the reduced graph is p — q, where 
g + 1 is exactly the incident number of the vertex w.r.t. the edges of the reduced graph. 
An example of such factorized tree T for p = 3 and g = 1 is presented in Fig. 1. 

Remark 2 Note that a set of proper closed paths of the graph T has the natural analogue 
in the continuous projective geometry case — it is a set of closed geodesies on the (open) 
Riemann surfaces. 

For each point ?/ G T, we define its distance to the reduced graph d{y, F) by the formula 
ini^^F d{x,y), which we sometimes abridge to d{y). This minimum is between the point 
y and a unique point x G -F, which we call an image t{y) G -F of the point y. Each branch 
B can be therefore naturally projected into its summit (in other terms, root, image) 
t{y) G F. 




Fig. 1. An example of the factorized tree T for g = 1, p = 3. 

Note that the terminology of ||2^, where the scattering theory on graphs has been 
explored from the Schrodinger potential standpoint, differs slightly from our notation. 
The reduced graph is the basis graph in Novikov's notation. Moreover, instead of infinite 
branches growing from vertices of F, half-infinite tails growing from the corresponding 
points (nests) have been considered. However, this setting is very close to the one used 
here because we consider (automorphic) functions that depend only on distance to the 
reduced graph in what follows; in this respect, points of tails just label these distances 
and the both approaches are equivalent in this respect .0 

2 Functions, operators, and potentials on graphs 



^An important difference is due to different types of the potential operators considered. 
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2.1 Automorphic functions and potentials on T. 

In the original setting by Bass [0], a unitary representation p : i^iiT) — > GL{vp) of 
the group V was fixed and a hnear space Vp was associated with any vertex (and, in 
principle, with any oriented edge) of the tree D. One can consider the space L{x,D,Vp) 
of automorphic functions on the tree D. Those are functions G{x) defined on V{D) 
(or L{D)), taking values in Vp, and satisfying the condition G(7x) = x{l)G{x) for all 
7 G r and x G -D. Here xil) is a left-invariant character of the group F. If F C -D is 
a fundamental domain of F in D then each G G L{x, D,Vp) is completely defined by its 
values on Y. 

However, we adopt an equivalent description, which is closer to physically meaningful 
lattice gauge theories. Namely, instead of considering automorphic functions, we consider 
periodic functions on D{F), which can be lifted to the graph T (the factorized tree), while 
nontrivial character set of the group F is ensured by introducing (unitary) potentials on 
the edges of the graph T. 

Let ^ be a group (or an algebra) and Vp the finite-dimensional representation of A. 

In 1^, the case of the trivial representation xil) = 1 was considered; nevertheless, 
already this simplest example manifests all the main features of the general theory. In 
the present paper, we demonstrate how the analogous method can be applied to the case 
of arbitrary unitary (in general, non-Abelian) representation {vp,x)- 

Definition 4 Given the group A and its finite-dimensional representation Vp, we associate 
the representation space v of the group A to each vertex and to each oriented edge of the 
graph T. Then, the corresponding linear spaces of functions defined on vertices or edges 
of a graph and taking values in Vp are Co = v^'^^ and Gi = f 1^' (the power is understood 
as the tensor product power). 

In what follows, we interpret tilde quantities as pertaining to the reduced graph or to 
its universal covering. For example, the linear spaces Gq and Ci correspond to vertices 
and edges of the reduced graph F. Note that the graph is always assumed to be dosed as 
the simplicial complex, i.e., every edge enters this graph together with its both endpoints 
(vertices) . 

We define also the natural bilinear forms on Gq and Gi: for f,g & Go, 

if, 9) = T. fixYgix), f{x),g{x) G Vp, (2.1) 

x&V 

and analogously for f,g G Gi. 

Definition 5 1. First type of potentials are the unitary potential on edges, Ug G A, 
where ^ is a (non-Abelian) group, the unitarity condition implies that Ug= UZ^ . Those 
are the potentials one encounters in the Bass case 
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2. Second type of potential are Hermitian potentials on edges lying in an algebra 
A such that A^ = A-^ = At (the reality condition) . 

3. Third, we call the nontrivial torsion potentials the potentials Upn) determined for 
the pairs of consecutive edges (/i precedes U) such that 

[/.I = U^^. (2.2) 

Another natural demand, which is often imposed on the set of torsion matrices Uftp 
is as follows. Let /ij, i = 1, . . . , n, be edges coming (in the fixed cyclic ordering) into a 
vertex of the graph. Then, for all 1 < i, j, k < n, 

U^^U^^ =U^^ ; (2.3) 
MiMfc MfcMj M»Mj 

these relations must hold for any vertex and they imply, in particular, that U^^ = h for 
any edge e G L{D). 



2.2 Hecke operators on graphs 

In this section, Co and Ci are the corresponding spaces of functions for arbitrary graph. 

We consider first the operators that act inside linear spaces Cq and Ci. All these 
operators are assumed to act on the corresponding universal covering D (or Dp) with the 
natural identification of (oriented) edges and vertices of D with their preimages in T (or 
in F). We denote by Jq and Ii the identity operators in the respective spaces Cq and Ci. 

1. We introduce the set of the Hecke operators M„ : Cq Cq, n = 1, . . . , oo, such 
that their action on basis vectors is 

Mnv,= Yl Ufl,---Ufl^Vy, (2.4) 

y- \Px,y\=n 

where the product runs over all oriented edges entering the path P^^y and x,y & D. 

For the tilde operators pertaining to simplicial complexes associated with the reduced 
graph F, the corresponding to ( |2.4D definition is 

MnV^= Ufl,---U^^Vy. (2.5) 

The action of operators ( |2.4| ) and ( |2.5| ) (as well as of all other operators) must be 
continued to the whole space Co by the linearity property. 

2. Next, we have the valency-counting operator Q : Cq ^ Cq, 

Qvx = (#neighbors - l)vx. (2.6) 

For the uniform tree, the operator Q is merely the identity times p. However, for a 
nonuniform tree, say, for Dp, the operator Q possesses some nontrivial properties. 
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3. The A-operators A('u) : Co — > Co and A('u) : Co — ^ Cq, 

A{u) = Iq-uMi+u^Q, A{u) = Iq-uMi + u^Q, m G C, (2.7) 

which have inverse operators (for \u\ < 1/p in the trivial character case), play an important 
role in the Bass construction. We sometimes label their representation dependence writing 
them as, say, Ap{u). 

Operators M„ ( p.4|) constitute a basis in the center of endomorphism algebra. Their 
multiplication algebra in the tree D is 

M,Mn = M^+,+pM^_i, n>2 

Ml Ml = M2 + (p+l)Mo, Mo = Io. (2.8) 
If if) is an eigenvector of Mi with the eigenvalue t, 

Miip = tip. (2.9) 
then it is also an eigenvector of all M„ and 

Mnip = Snip,t)lP (2.10) 

where Sn{p,t) is a system of orthogonal polynomials in t with a generating relation 

tSn{p, t) = Sn+l{p, t) + pSn-l{p, t) , Si = t, S2 = + p + I. 

In the general case of a nonuniform tree Dp, the operator Q obviously does not 
commute with the Hecke operators, which in turn also become mutually noncommutative. 
Nevertheless, it is possible to write algebraic relations analogous to ( p^.Sj ) using the operator 
Q: 

MiM„ = M„+i + gM„_i, n>2 

MiMi = M2 + (Q + l)Mo, Mo = id. (2.11) 

Moreover, 

oo 

5] «"M„A(u) = (1 - n2)/o. (2.12) 

71=1 

Relation ( p^.l2D implies that the operator A(u) is the generating function for the operators 
M„. 

4. The Laplacian of a unitary theory is merely A(l). For this Laplacian to be a 
Hermitian operator w.r.t. brackets ( |2.1|) we must impose the unitarity condition on the 
potential, 

Ut = U^\ (2.13) 

In the Laplacian was determined as the operator of the Schrodinger type, i.e., 
it has form ( |2.7D with the first Hecke operator Mi as in ( |2.4| ) but with the Hermitian 
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potential Ari = (elements of an algebra) substituted for Urt (elements of a group) on 
the edges of a graph. The corresponding Laplacian 

y-\P^,y\='^ 



is also Hermitian, but no relations of type ( p.ll ) exist 



5. On the space Ci we first define the inversion map J : Ci ^ Ci, which merely 
changes all orientations of edges, 

Je = Use, (2.15) 

Moreover, there exist the set of Hecke operators for the space Ci. They all are generated 
by a single operator T = Ti : Ci ^ Ci. 

Definition 6 The first Hecke operator acting on the space Ci is 

Te= ^e-e^. (2.16) 

following e 

Then, we can obviously define T„ as the sum over all terminal edges of oriented reduced 
paths of the length n + 1 starting from the given edge, 

T'^{eo)= E (2.17) 

(eo,ei,---,e,„)red 

However, in contrast to the space Co, the relation between T„ and Ti is merely T„ = TJ^ 
for any graph, and the family of these operators is commutative for every tree, no matter 
uniform or nonuniform. Then, for any graph, the corresponding generating function is 
1 1 - uTi. 

CO 

E m"T„(/i - uTi) = Ji, ueC. (2.18) 

n=0 

Definition ^ can be easily generalized to other potentials from Definition ^. For f/^; 
replaced by the Hermitian matrices Af^, we just obtain the Novikov ||26| potential on edges. 
Important examples are provided by potentials of the third type (with the nontrivial 
torsion). Then, in the most general case, we can define T as 

Te= E UgFgg.e^, (2.19) 

Cfj, following e 

where neither Ug nor Fgg^ must satisfy the unitarity or Hermiticity conditions; in formula 
( p.l9| ), the potentials on edges with opposite orientations can be set arbitrary, the L- 



function expression through the operator T, which are presented in the following section, 
remain nevertheless valid. 



We note two important examples of ( ^.19| ). The graph description [^, of the 
Teichmiiller spaces of complex Riemann surfaces discussed in Sec. 6; the other is the Ising 
model on a lattice governed by the Kac-Ward operator determinant |^. 



14 



2.3 Intertwining operators between Cq and Ci. 



Following 0, we introduce a set of operators acting between two spaces Co Ci. 

1. First, we have two boundary mappings di,dQ : Ci Cq, 

die = AsXi (2.20) 

and 

doe = xo, (2.21) 

2. We also have the coboundary operator ctq : Cq ^ Ci, 

\Px.y\ = l 

which set into the correspondence with the point x a linear combination of all edges 
outgoing from this vertex. (In 0, the operator (Xi, which sets into the correspondence to 
the vertex x all incoming edges was defined. However, this definition, which is meaningful 
for a nonpotential case, becomes ill defined when a potential is introduced. Fortunately, 
in what follows, we do not need the operator cti.) 



2.4 Relations between operators 

Let us introduce the total (half-) dimensions of the spaces Cq and Ci, 

ro = rank (Cq), = rank (Ci)/2. (2.23) 

We also introduce auxiliary operators 

d{u) = DqU — di, a{u) = a^u. (2.24) 

The standard boundary operator is then d ^ — di = 9(1). 

Direct calculations show that the above operators satisfy the following relations |@|, 
which hold for any graph: 

doCFo = g + 1; 
^i^o = Mi; 

T = cTo^i - J] (2.25) 
d{u)a{u) = A{u) - (1 - u'^)Io; 
a{u)d{u) = u{T + J){uJ - h). 

An important assertion proved by Bass establish a connection between the spectral 
Laplacian problem and the spectral problem for the Hecke operator Ti. 
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Lemma 1 [0 For the potentials of the first type in Definition we obtain 



det(/i — uTi, 



Co 



Proof. Let us define the matrices acting in Co © Ci, 

,2^ 







Then from ( ^.25 ), we obtain 



LM 



A(n) 
(r{u) 







d{u) 

h 



and ML 



M 



-d{u) 
'l-u^)h 



a{u){l 



u 





uT){h 



uJ) 



The assertion of Lemma |l| follows from equating the traces of LM and ML; we only need 
to evaluate the determinant of Ji — uJ . It is block-diagonal in the basis of (nonoriented) 
edges. Each edge admits two orientations, so we have 



det(/i — uJ) = det 



h 



-uU 
h 



^ \En\ 



(1 



u 



2\ri 



[2.2Q) 



for any matrix U . Note that ri is exactly the dimension of the representation of the 
group U times where \E\ is the number of unoriented edges. Therefore, tq — ri = 
dim ?7 X — lE'l) where the difference between the numbers of vertices and (nonoriented) 
edges is exactly 1 — g — the Euler characteristic of the graph. 

In general, no relation of type of Lemma |I| exists for potentials of the second and third 
type from Definition ^. Therefore, the spectral problems for the operators A(m) and T{u) 
are not related in these cases. However, for graphs of the third kind from Definition ^ with 
conditions ( |2.2| ) and ( p.3| ) imposed, we can, nevertheless, connect these two operators for 
the price of "blowing up" the vertices of a graph; this construction is discussed in Sec. 6. 



3 L-function 

Definition 7 [Q] Let Vp be a C'^-module with the character xp • T{C\T]) C^. Then, 
the IharaSelberg L-function L{p,u) is 

L{p,u)= n det(/,-M^("V(^))-\ (3.1) 

where the product ranges all conjugate classes of primitive elements of the group F, or, 
equivalently, all proper oriented closed paths without periods in the graph T, are lengths 
of these paths, and p{^w) = Ug-^ . . . Ug^^ is the product over the path that corresponds to 
the element G F of elements of the group A. 



We present now the Bass' formulation and proof 0] of the Hashimoto's theorem . 
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Theorem 1 Let p : F — * GL{V) be a finite- dimensional complex representation of T. 
Then, the L-function ( |3.1| ) is a rational function in the variable u: 



L{p,u) 



1 



det(/i — uTi), 



(3.2) 



where Ti is from ( p.l6| ) . 



From Lemma |1] and Theorem |1] the corollary follows 



Corollary 1 // the unitarity condition ( |2.13| ) is satisfied, then the L-function can be 
expressed through the determinant of the operator Ap{u): 



L{p,u) 



-1 



2\ro-ri 



det(Ap(M)). 



(3.3) 



The product in ( |3.1D converges absolutely in the domain {u : |m| < (gma^l-^l)^^}, 
where Qmax < p is the maximum incident number of the tree Dp and |^| is an absolute 
value of a maximum element U-^j. Also, we have for the logarithmic derivative of L{p,u), 



where the sum ranges all conjugate classes {7} distinct from {1}. 
Proof of Theorem 1. 

Let us choose in the tree D some fundamental domain 'P(F) of the symmetry group 
F. We now consider the contribution to the trace of the operator T"* coming from some 
edge e*! G ^^(r)- Only those e j G -D contribute to trT"* that 

1. lie on the distance m to the initial edge Ci (along some reduced path 6162 . . . e^+i); 

2. have the orientation along this path. 

The crucial observation follows. Let us consider the action of some element 7 G F on 
oriented edges of D (Fig. 2). It is easy to see that if the edge does not belong to the 
invariant axis D{'y) of this element (like Cx does not on Fig. 2), then the orientation of its 
image under the action of the element 7, 762;, is opposite to the orientation of the reduced 
path coming from to 'ye^- Therefore, no such edges contribute to the trace! Only the 
edges that (like Ci in Fig. 2) belong to the invariant axis .0(7) preserve their orientations 
toward -D(7) under the action of 7 (the translation along -D(7)). Say, for ci, we have 
ei+i = 761 in Fig. 2. 



d 

du 



log(L(p,«)) = 5:A(7)tr(p(7))n'W-\ 



(3.4) 
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ei ei+i=7ei ^(7) 

Fig. 2. The action of tlie group element 7 G F on tlie edges of 
witli tlie length /(7) = /. Note that 76^: 7^ T''""*"^ei for all m. 

Note. This is why we use trT™ instead of, say, Mm — the mth Hecke operator for 
the space Co — for trM^, all pairs of vertices x,y & D, not only those belonging to the 
invariant axis -D(7), contribute to this trace as soon as d{x,y) = m and y = ^x. 

Coming back to the space Ci, we fix for a moment the edge e\ G -D(r) C D. As was 
mentioned above, for each primitive conjugate class w one can find not necessarily one 
axis D{w) C Dp such that Ci G D{vj). We denote by the oriented edges of the reduced 
graph F itself, thus each Cj G -D^ is an image of the edge & F with the orientation 
naturally preserved. For each element 7 G F, we then set into the correspondence the 
periodic sequence of "letters" (the cyclic word) 

£)(7) ^ . . . {dia2 . . . di){aid2 ... a;) ... , (3.5) 

where I = A (7) = l{w) is the length of the generating element for 7. 

Note that the elements of F determine the symbolic dynamics of "words" [Lyndon 
words, see [|l^) composed from the "letters" aj, which can be noncommutative elements 
of the group A. We even do not need a unitarity condition ( |2.13|) in this case; in particular. 



Theorem |l| holds even if letters a, and Oj are not related. Thus, for the rest of the proof, 
we merely identify the letters di with the elements Ug-Ci and the proof is valid for any 
(non-Abelian) group or algebra elements dwelling on the oriented edges of the graph F. 

It is clear that if ^(7) = m, then l{w) = I, l\m, that is, I is a divisor of m. Note 
again that there can be repeated di in the cyclic word 0102 .. . di{di+i = di), Moreover, 
this sequence may contain subperiods if 7 is not a primitive element, but, obviously, the 
minimum length of this subperiod, /, is a divisor, l\m. 

We now fix the element di — the preimage of ei. In order to determine its contribution 
to tr T"^ we must find all possible different cyclic expressions 

did2 . . . dm{dm+i = di) (3.6) 

that include this element. We now establish a correspondence between sets (|3.5|) and 

(PD- 
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1. To each finite cyclic sequence ai . . .0^(01) we inambiguously set into the corre- 
spondence the infinite periodic sequence . . . (ai . . . am){ai . . . am) ■ ■ ■ corresponding to a 
unique element of primitive conjugate class {w} with l{w) = l\m. 

2. On the contrary, we now choose an element from {tJj}, or, equivalently, from some 
periodic reduced sequence 

...(.. .){ai^ai2 . . . di^){ai^ ...)... (3.7) 

with no subperiods and with the minimal period / = l{w). If it contains the (oriented) 
edge 0,1 di times among the edges {a^^, . . . ,d^} and, moreover, / is the divisor of m, then 
there are exactly di different sequences di . . . dm{di) (|3.6|) containing in ( p.7| ). Eventually, 
doing the sum over all edges of F (or of the fundamental domain 'D{T) G -D), i.e., 
evaluating the trace of the operator T™, we find that the contribution from the sequence 
...(.. .)(aij . . . J . . . to this trace is 

\Lf\ 

^ # {aj in {ttii, . . . ,aij} = /, (3.8) 
i=i 

where / is the total length of the primitive element w. Therefore, we have the following 
remarkable formula: 

iTT^ = Y.i- ^Ap{^r% (3.9) 

l\m {ro:i(ro)=/} 

where the sum runs over all primitive conjugate classes of F. 

We now obtain from (^.4|) the representation for the L-function: 



J 00 
u-\og[L[p,u)) = 5:5:/(t^)-«"'(-)tr[p(^)T 

{u7} n=l 

°° d 
= from(3.9) ^ trT^-u'' = -M— logdet(l-u-T). (3.10) 

k=i 

This completes the proof of the theorem. 

4 Spectral problem on infinite graphs 

We set the spectral problem considering a function ip E Cq such that 

L-^ = Xip, A e C, (4.1) 

where the Laplacian L is 

L = M,-{Q + 1).. (4.2) 

We first determine the solution of ([4.1|) on a branch (tail). We assume the graph T to be 
uniform {{p + l)-valent) but infinite. Then, we consider a branch (tail) of the graph T 
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with the starting vertex xq (the index i is the distance to xq in the reasoning below). We 
have ^'(a^i*^) = v^^^; then, 

E + -iP+ 1)^* = ^ > 0' (4-3) 

s=l 

where x\'}_i is a unique vertex (from the {i — l)th layer) that precedes Xi (when moving 
along the branch from its root). 

The general solution to ([4.3|) is 

^f^ = E«±(4^-rV^.-(4"yM, (4.4) 
± 

where the product ranges all edges of the path P (3.), f° and f° are constant vectors 
from V, and a± are the solutions of the equation 

pa|-(p + l + A)a± + l = 0; t = ^^^^^ , a± = t±^t'^- l/p. (4.5) 
The vectors 

[vf\], = «i(f^^r^ V' • • • {u'rWvD^, (4.6) 

where (f^)^ is the pth component of the vector remain collinear for any i and their 
ratio (a+/Q;_)*(f° )p/(f° )p depends only on the level i. Therefore, it is natural to relate 
the data {v^j^)p and {v^_)p with the scattering matrix data on the graph T. For this, we 
call solution ([4.4|) with only v° nonzero the incoming wave and solution (|4.4|) with only 
v° nonzero — the outgoing wave. 

Here it becomes clear why we can relate the branch and tail descriptions. First, we 
can perform a gauge transformation in order to eliminate the dependence on potentials 
on edges of branches. For this, we merely set 

No dependence on potentials on edges remains then in the transformed variables, and we 
obtain merely that 

N±]p = a±bo±]p 

for any path of length i starting from the summit of the branch. In what follows, we set 
Upi = Iv on all external {i.e., not entering the minimum reduced subgraph) edges. 

We can now identify all vertices and edges that lie at the same distance from the 
summit setting 

vf p^'-^^'^v, for s > (4.7) 

and erasing the now redundant superscript {i) (see Fig. 3). We thus obtain just the 
picture of where a finite (probably zero) number of tails (half-axes) begin at a vertex 
of the basis (or, in our notation, reduced) finite graph. We make transition ( [4.7|) in 
order to preserve the scalar product (|2.1| ). Simultaneously, we must scale and shift the 
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corresponding eigenvalue, A — >• [A + {^/p — only after this operation, the answers 

for two branches of solution ( [4 .51 ) in our approach and in the approach of do coincide. 




^^0 
Vl 



Fig. 3. The branch — > tail transition on the external edges of T. 



Definition 8 Let the boundary points Xi G OF are those points of F that have nonzero 
number of external tails starting at these points. We also consider the natural embeddings 
of linear spaces ((?F)o,i C Co,i (recall that (-F)o,i = C'o.i)- 

4.1 Spherical functions on factorized trees 

We now introduce spherical functions for the graph T. First, we consider them on the 
tree D itself. Choose the vertex x G V{D) and claim it the center of the tree. Then 
the spherical function ip{n, x) is an eigenvector of Mi, Miip{n, x) = tip{n, x) that depends 
only on the distance in the tree from the point x, i.e., it is constant on each sphere 
S{n,x) = {y E V{D) : \Px,y\ = n}. In the case of trivial potential (trivial representation 
of the group F), we have 

i>{n, x) = a+a" - a_a" , (4.8) 

where (cf. (^)) 

= 1/p, 

and setting ?/'(0,x) = 1, we obtain 

pa+ — a- pa^ — 



(p + l)(a+ - a-) (p + l)(a;+ - a;_) 



(4.9) 



In order to define a similar object for a general graph T whose "center" now is the 
reduced graph F, we consider a superposition of solutions (4^) with the sources Sy placed 
at the vertices of Dp, 

^jj(x) = Sy^P{d{x,y),y), (4.10) 

such that the function Sy is periodic w.r.t. the action of F: s^y = Sy for all 7 G F and 
y E Dp. The behavior of this solution on each branch "growing" from dF is described by 
( [4. 81) with the preexponential factors depending only on the point of dF into which this 
branch can be projected. 
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Let us introduce a kern function K: 



K{z,y\x) = Y.x''^^^^^y^^ X n Ufl^^ (4-11) 

where z, y are points of Dp. This function is periodic under the action of F over both its 
arguments separately. Therefore, it is well-defined on T itself. For the trivial representa- 



tion X = 1, this function is also symmetric in z and y. Function ( [4. 101) then becomes 



iIj{x) = ^ Sy[aj^K{t{x),y\a+)a'il:^^ — a_K {t{x) , y\a^)a'5^^], (4-12) 
yeF 

where t(x) G -F is the image of the point x E D. (If x G -F, then t{x) = x.) We call the 
part of ( |4.12| ) proportional to a'i^^^ the retarded wave function and the part proportional 
to the advanced wave function. Therefore, tplx) has a general form 

^{x) = A,dMx))af ^ - = - (4.13) 



Proposition 1 Given a set of vectors v^{xi), Xi G OF, the set of v^_{xi) is uniquely 
determined for allt E C except a finite set of points {the points of the discrete spectrum). 



The important case of functions (|4.13| ) is where ip-\-{x) is nonzero on only one branch 



growing from a single point of dF. Then, we can define the scattering matrix as follows. 

Definition 9 The scattering matrix (S-matrix) S{t) is the square matrix of the size 
rank (9F)o with the entries Si^j^, \{ia}\ = \{ji3}\ = rank(c}F)o, Si^j^ is the value of 
(f° )^(xj) for (f° )^(xfe) = Sa'ySki and {v'^) pertain to the solution of spectral problem (^4.1|) 
on the whole graph, i.e., entries of the S'-matrix correspond to solutions of ( |4.1| ) of the 
form 

{jp} 

We postpone the detailed study of the spectral properties of the S'-matrix to the next 
section and formulate here the main theorem connecting S'-matrix and determinants of 
the operators on the reduced graph. 

Theorem 2 In the case of unitary potentials on edges, 



det S{t) = 



^^yank(F)o det(A,(a_)) 



I det(Ap(a;+)) 



where (Ap(a_|_)) is operator (|2.7] ) in which we explicitly indicate the dependence on the 
representation p. 
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Proof. Each solution of spectral problem ( [4.1|) having form ( |4.4| ) outside the reduced 



graph can be split into the advanced and retarded partial wave functions in the total 
graph (see ( |4.13| )) where 

^±(^)= E E «±"±^"'^'^ n U^'-<^o.). (4.15) 



Here s{xa) : Fq ^ V is the source function and a± are from (|4.9|) . 

Note that the "splitting" into the advanced and retarded parts or, equivalently, the 
source function s{xa) is inambiguously defined everywhere outside F (the splitting is also 
inambiguously defined at boundary points x G OF), but is ambiguous at the inner points 
of F. 

In the space of s{x) (in fact, this space is Cq), there exists a linear subspace of source 
functions that generate one and the same eigenfunction iIj{x). For example, if Xq G F, 
xo ^ dF, we can propose two source functions si{x) = {0,x ^ xq\ \,x = xq} and 
s^ix) = {0,d{x,xo) 7^ 1; 2p/t,d{x,xo) = 1}. Then, these two source functions generate 
the same eigenfunction iIj{x). Moreover, for almost all A and s{x) G Cq, there exists a 
unique s{x) G Co such that both s{x) and s{x) generate the same ip{x) and s{x) = in 
all inner points of F. 

Therefore, the problem is to define a "convenient" splitting ( [4 .131) (an analogue of a 
finite-dimensional "gauge fixing" in field theory). 

Let us consider eigenvalues Cn of S(t). Those are numbers for which there exist solu- 
tions of ( [4.1|) with the sets (f°)a(xj) and (f°)a(a;j) proportional to each other: 

{v^)a{xi) = c„(f° )a(xj) for all a and Xj G dF. (4.16) 

We now expand the condition ( |4.16| ) to all (external as well as internal) points of F. 



Definition 10 Let the spherical function on the graph T be the solution of spectral 
problem ([4.1|) of form ([4.15|) such that 



[^+(:^)U[^-(^)L = c„(^) (4.17) 

where Cn is the constant independent on the point x and the representation index a and 
dist [x, F) is the well-defined distance between the point x and the reduced graph F (for 
X G F, dist {x,F) = 0). 



Lemma 2 



det S{t) =Y[c„ 

n 



Actually, condition ([4.17|) fixes the choice of the admitted source functions s{xa). 
At the same time, this is a system of linear homogeneous equations, which has nonzero 
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solution only for a finite set of c„. Below, we present this system in terms of the Hecke 
operators ( |2.5| ) on the graph F. 

Considering Eq. ( [4.17| ) at the points of the reduced graph, we obtain 

or 

^ Sya+K{z, y\a+) = CnYl Sya_K{z, y\a_). (4-18) 

Then, obviously, 

„ detK{x,y\a+) /g+V'^""^^^" 

where K{x,y\a) is the kernel for the linear operator acting on Co (note that because F 
contains all loops of the graph T and is connected, then, if x,y G F, then any path 
Px,y C F): 

(^^)(^) = EE«"'^'^' ( n uAs{y). 
Considering the universal covering Dp, it is easy to verify that 

K{a)K{x,y\a) = (1 - a^)/^ 

and, therefore, 

det K{x, y\a) = (1 - a^Y° det -^A(a). 

Now, using the identity 

a+(l — a^) a+ 
a_(l — a?^) q;_ 
and formula ( [4.19| ), we obtain the assertion of Theorem |^. 

Theorems ^ and ^ imply the following correspondence between the scattering matrix 
determinant and the L-function of the reduced graph: 

detS(i)=f^)7^)'"'^^, a.a- = l/p. (4.20) 

where L(p, a) depends only on characters of loops of the reduced graph F. Another form 
of ( F20D is 

det5(t)=f^y°^^^5^. (4.21) 
\a-J detA(«+) ^ ' 

Note that only the first "volume factor" in ( [4.20| ) and ( [4 .211) depends on the volume 
of the graph F, all nontrivial factors depend only on loop characteristics of F. 

Since = 1/p, the combination {a±Y° det A(Q;:p) acquires the form 

{a±Y° det A(a^) = det(a±A((a^)) = det la± + - ^mA . (4.22) 
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Remark 3 Returning to ( [4.18|) , we see that the determinant of the operator K{z,y\a±) 
is nondegenerate (everywhere except discrete points of spectrum of A(m)). It is in contrast 
to the aheady mentioned ambiguity in choice of the source function s{x). Let us now act 
by the operator A(a+) on both sides of ( [4. 181 ). For any interior point of F, 

(a_A(a+)s)(x) = (Q;+A(a;_)s)(x); 

therefore, we obtain siy) = Cis{y), i.e., s{y) = in all internal points for all q ^ 1. 
Hence, using condition ( [4.18| ) we gauge out the whole ambiguity due to the choice of s{y) 
at almost all values of the spectral parameter t. (The cases where such a condition fails 
just correspond to the exceptional points of the discrete spectrum, see below.) 



4.2 Structure of spectrum 

We now study the spectral properties of problem (|4.1| ). These properties turn out to be 
quite similar to the spectral properties of Schrodinger type potentials investigated in fST 

We now consider eigenfunction problem ([4.1|) as the spectral problem. 



Definition 11 The continuous spectrum, or the scattering zone of the spectrum, is the 
domain < 1/p where a+ and a_ are complex conjugate to each other and their absolute 
value is exactly l/i/p; then, 

A G (-P - 1 - 2^^, -p - 1 + 2^). (4.23) 



The points of normal discrete spectrum are eigenvalues of problem ( [4.1|) that lie out 



side the continuous spectrum domain ( |4.23| ) and correspond to real solutions that decay 
exponentially (as a^) with |a+| < l/y^ and such that these eigenfunctions are nonzero 
at least in one branch. 

The points of the exceptional discrete spectrum may appear for any real A. Those are 
points such that the corresponding eigenfunctions vanish identically outside the reduced 
graph. 

The first two types of points are customary for a majority of spectral problems. The 
third type is rather specific for graph problems or, more generally, for discrete spaces (its 
existence in the problem under consideration was observed in pO[|). 

Note that points of the exceptional discrete spectrum exist not for all graphs; cus- 
tomarily, it is some graph symmetry (e.g., the Z2 symmetry, see |2^) that is responsible 



for the appearance of such points. Obviously, for such an exceptional solution to exist, it 
must vanish not only on all branches, but also at all summits of these branches, i.e., at all 
boundary points x G dF (and, correspondingly, at all boundary points of the minimum 
reduced graph). This condition is not easy to formalize; however, in the unitary potential 
case, we can formulate some regular criterion for the existence of exceptional discrete 
points in spectrum (see Proposition |^) . 
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If k is the total number of branches growing from the points of OF, then, for any 
A from the scattering zone, considering the spherical function case (the functions that 
are constant on spherical slices of branches) we always have a /c-dimensional subspace of 
eigenfunctions corresponding to the given eigenvalue A. 

Fixing the spectral parameter A, we can introduce the 2A;-dimensional space H'^'^ of 
asymptotic states — the asymptotic spherical function-like local solutions of ( |4.1| ), which 
always has two independent solutions proportional to a" and a" in each branch. These 
functions, ijjf, i = 1, . . . , k, constitute a basis in H^''.. Then, the true scattering problem 
solutions (|4.14| ), i.e., those that can be continued to the whole graph, span a /c- dimensional 



2k 



subspace A'^ of H 

Problem 1 Is it possible to develop a Lagrangian plane description for the unitary po- 
tential case in analogy with the Schrodinger potential case [P^]? 



Proposition 2 The S -matrix ( |4.14| ) is unitary in the scattering domain ( [4.23[ ). 



Proof. The proof resembles (however, with some important variations) the proof of 
the analogous statement for the Schrodinger potential scattering case. 

As the first step, we define the symmetric function W{(j),ip) (we call this function 
s-Wronskian by analogy with the antisymmetric Wronskian in [Q), which sets into the 
correspondence to two functions from Cq the function from Ci by the formula 

W^(0,^)= i<P*i^)Ufi.,y^{y) + rix)U^.,ym)- (4-24) 

If and ip are two eigenfunctions from A'^ (by definition corresponding to the same value 
of the spectral parameter A), then the s- Wronskian is a chain (one-cycle with possible 
open ends), i.e., in the infinite graph, we have 

dW{(j),ilj) = 0. 

This follows from the simple calculation, 

51^(0, ^)(x) = 4,*(x)J2U^^,Jiy)+rix)J2U^^Jiy)- 

y y 

y y 
= 0*(x)(MiV^)(x) - (0*(x)(MiV^)(x))* + 

+ilj*{x){Mi(t)){x) - (^*(x)(Mi0)(x))*. (4.25) 

Using now the relation Mi0 = (Q + 1 — A)0 and remembering that Q and A are Hermitian, 
we obtain zero in the r.h.s. of (|4.25| ). 

The s- Wronskian of two solutions is therefore a chain. This means in turn that if we 
now cut the branches of the (infinite) graph at some distance n from the reduced graph F 
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(this distance must not be the same for all branches; it is only preferable to be the same 
on each branch), thus obtaining some finite graph Tcut, then, nevertheless, the total sum 
vanishes: 

J2 dW = 0. (4.26) 
Here the sum ranges all boundary points x G dT^^^, and dW is restricted to the graph 

It is convenient to introduce the bilinear pairing on H'^'', 

on the boundary of Tcut for two arbitrary asymptotic vectors from H^'^. 

The reasoning in the beginning of Sec. 4 show that we can consistently eliminate the 
gauge potential dependence on all external branches. Then, the ipf functions are mere 
exponents of a±: tpf = a!^ on ith branch and zero otherwise. Then, using the identities 
= 1/p and a_ = which hold in the scattering zone, it is easy to find that 



2p 

and 



Now we remember that for two arbitrary solutions ip and of spectral problem ( [4.1|) de- 
termined on the whole graph, the expression for dWlip, (p) exactly coincides with {{ip, 0)), 
where only asymptotic states must be taken into account. Consider now the eigenfunc- 
tions = ip^ — J2j Sijip^ (cf. formula ( [4.14|) ).p| Then, for two such functions, we 
have 

= ((VI/,, vl/^.)) ^ ((^+, ^|J+)) (6,, - J2 . (4.27) 

Because {{ip^ ^ip^)) ^ 0, we obtain that the ^-matrix is unitary. The proposition is 
therefore proved.0 

It is interesting that we can deduce the existence of the exceptional discrete spectrum 
from the L-function itself. 

Proposition 3 Exceptional discrete spectrum appears iff^a+,a- : a+a- = 1/p and the 
L-function has poles at both and a_. 



Sketch of the proof. We seek the exponential solution of spectral problem (|4.1| ) such 
that one of the waves, retarded or advanced, just vanish identically in all vertices of the 



■^We use condensed multiindex notation. 

''Considering formula (4.20) we obtain that the 5- matrix determinant is a unitary function as soon as 



a± he in the scattering zone; this is of course in accordance with Proposition |^. 
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graph T (including now also internal points of F). This immediately implies that there 
exists a vector (source function) /o G Co such that 

A(a+)/o = 0. 

Then the nonlocal operator K, 

Kfix) = J2Kix,y\a)f{y), f e Co, 
yeF 

which is inverse to A (a), becomes ill defined. 

Let fo be the null vector of A (a). This immediately produces the characteristic 
equation on a: 

det A(a+) = 0, 

and if simultaneously A(q;_) 7^ 0, then, following Theorem H, the S'-matrix develops a 
singularity (zero or pole), which indicate that for such value of A there exists a solution of 
spectral problem ( [4.1|) in which one (and only one) of the retarded and advanced waves 
is nonzero (cf. |]30|). 



Thus, if no cancellations between numerator and denominator occur in formula ( 4.20 ), 
i.e., the assertion of the proposition is not true, then the discrete spectrum is exhausted 
by regular eigenfunctions and no exceptional spectrum exists. 

On the contrary, when there are such poles of the L-function that the assertion of the 
proposition is satisfied, then the regular discrete spectrum does not provide all possible 
solutions of the characteristic equation; the solutions that yield zero both in the denomi- 
nator and in the numerator of ( [4.20| ) must therefore correspond to the exceptional discrete 



eigenvalues. 

Problem 2 To make this sketch of the proof the formal proof. 



5 Examples. Algorithm for calculating L{l,u). 

For simplicity, we assume in what follows that F = D(r)/r. We omit most of the proofs 
that are simple exercises in linear algebra. 



Example 1 Graph T with no external branches. Here we consider a limiting case of 
our construction when the graph T exactly coincides with F. In this no-scattering case, 
q[z) = p and we obtain the inversion relation for A(m). If a.a^ = 1/p, we have 



det"^[(l+pa+)/-a+Mi] = det 



1 + 



1 



1 



-Ml 



pat pa. 
alpY' det ~i[(l + pal)! - a_Mi] 



a 



ro 



— det "^[(1 + pal)! - a^Mi] 



a 



(5.1) 
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and using ([4.21|) and ( ^.1|) , we obtain that the total scattering matrix C = 1. Therefore, 
C is trivial and does not depend at all on the shape of the reduced (or, in this case, total) 
graph T = F. This is in accordance with the fact that all spectral points are points of 
the exceptional discrete spectrum in this case. 



Example 2 A one-loop case {g=l). In the case where the graph T contains a single loop, 
the L-function in the case of the trivial character xi"^) = 1 is 

^(l'^W= (T^^' (5-2) 

because the group Fi contains only two primitive elements 7 and 7"^ with the same length 
/(7) = /(7-') = n. 



5.1 Two-loop case. The general structure of L{l,u). 

We present an algorithm for calculating L{u) = L{l,u) for arbitrary graph. Let k be 
the lengths of edges of some graph F and = Then, from Theorem 1, we have 
L^^{u) = det(l — uTi) = P2ri — the polynomial of degree 2J2h in the variable u. This 
polynomial can be obviously treated as a polynomial in (all closed geodesies are built 
from intervals connecting three- and more- valent vertices of F) and we may treat 
as formal independent variables. Assuming that one of h is greater than the sum of all 
others, we conclude that the maximum degree of the polynomial P2ri {u) in each Oj is less 
or equal 2 Thus, L~^{u) = P2{{(ii}). This polynomial has integer coefficients and from 
the relation 

L-\u) = det(l - uTi) = (1 - u^f'^ det A{u) 

it follows that L~^{u) has zero of order at least (7 — 1 at the point u = —1 and zero of 
exact order g at u = 1. This follows from the relation det A(m)|^^-^ = 0, which holds for 
arbitrary graph, and from the relation 

d 

— det A(m)| -, = # maximum trees, 
ou 

(see Q where it was noted that the derivative of A(?i) at m = 1 is expressed by the 
Kirchhoff formula through the total number of maximum connected trees in the graph F 
and is therefore always nonzero). 

We turn now to two-loop case. There we have three possibilities (see Fig. 4, a, b, 
and c). 
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a b c 

Fig. 4. Three possible graphs for g = 2. 

The analysis in |Q shows how one can easily calculate the corresponding L-functions 
from the above considerations. The answer is 

case a = 1 ~ 51 ~ 2010203 1 ~ X] o.iO'j + 201O2O3 , 

\ i<j J \ i<j J 

case b = (1 — oi)(l — 02)[(1 — oi)(l — 02) — 403O1O2], 

case c L^J = (1 — oi)(l — 02)[(1 — Oi)(l — 02) — 401O2]. (5.3) 

(Note that the case c can be obtained either from a or from b case by setting 03 = 1.) 



Example 3 Comparing the Hermitian and unitary potentials. We now consider L-func- 
tions and the corresponding determinants for two graphs depicted in Fig. 5 endowed with 
Abelian potentials. 



^ — ( 


» — .^^^ 






1 — .^^^ 


Oi 

^ — ( 


» 


O3 ' 


V c 

— 1 




1 — ^ 


a 






b 





Fig. 5. Examples of reduced graphs with Hermitian and unitary potentials. 

We present answers for Hermitian and unitary potentials (where Oj — > o~^ for the 
inversely oriented edges) in case a. For the Hermitian potential, we obtain 

det A(m) = (1 + 2u^)2 - s^M^ s = a + b + c, 

det(l - T{u)) = 1 - 2u'^{ab + be + co) + u^{abb + be + caf - 4u^a%'^c^, 
while for the unitary potential, these quantities are, of course, related: 

det A{u) = (1 + 2u'^f - ssu'^, s = a + b + c, s = o"^ + 6"^ + c"\ 
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det(l -T(m)) = (1 -M2)det A(m). 

Case b in Fig. 5 is of interest because it is just the case where the exceptional discrete 
spectrum is nonempty for p = 2 and for any a (it suffices to set equal in the absolute 
value positive and negative charges at the upper and lower vertices in Fig. 5b; then, the 
eigenfunction vanishes identically at the middle point). For the L-function, we obtain 

L{u) = {1+UX + 2n^)(l -ux + u"^ - u^x + 2u^), x = a + a"\ p = 2, 

i.e., the condition of Proposition]^ is satisfied for any a, because the product of two roots 
of the first quadratic polynomial isl/2 = l/pin this case. 



6 Teichmiiller spaces via graphs 

In this section, we consider Teichmiiller spaces — the spaces of complex structures 
on (possibly open) Riemann surfaces S with holes (punctures) modulo diffeomorphisms 
homotopy equivalent to identity. In the vicinity of a boundary component, the complex 
structure can be isomorphic as a complex manifold either to an annulus (hole) or to a 
punctured disc (puncture). 

The graph description following and Jl^ is suitable for considering the finite 
covering T^{S) of the Teichmiiller space T^{S). A point of T^{S) is determined by a 
point of T^^S) and by the orientation of all holes of S that are not punctures. (This 
covering is obviously ramified over the subspace of surfaces with punctures.) 

It is well known that an oriented 2D surface with negative Euler characteristic can be 
continuously conformally transformed to the constant curvature surface. The Poincare 
uniformization theorem claims that any complex surface 5* of a constant negative curva- 
ture (equal —1 in what follows) is a quotient of the upper half-plane H+ endowed with the 
hyperbolic metric ds"^ = dzdz/{Qzy over the action of a discrete Fuchsian subgroup A(5') 
of the automorphism group PSL{2, R), 

S = H+/A(S). 

In the hyperbolic metric, geodesies are either half circles with endpoints at the real line 
R or vertical half-lines; all points of the boundary R are at infinite distance from each 
other and from any interior point. 

Any hyperbolic homotopy class of closed curves 7 contains a unique closed geodesic of 
the length /(7) = |logAi/A2|, where Ai and A2 are (different) eigenvalues of the element 
of PSL{2,'R) that corresponds to 7. 

Since Strebel ||3^, the fat, or ribbon, graphs have been used to coordinatize the Te- 
ichmiiller and moduli space. We use a rather explicit and simple version of this descrip- 
tion m . 
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Claim 1 For a given three-valent fat graph T of genus g and number of punctures n, 
there exists a one-to-one correspondence between the set of points ofT^{S) and the set 
edges ^fj^g^g of this graph supplied with real numbers (lengths). 



We propose the explicit way how to construct the Fuchsian group A(S') C PSL{2, R), 
which corresponds to a given set of numbers on edges of a graph T G T{S) such that 
S = H+/A(S').0 For this, we must associate an element G PS'L(2,R) to any element 
of the fundamental group 7 G 7ri(S'). 

To each edge a we associate the matrix G PSL{2, R) of the Mobius transformation 

In order to parameterize a path over edges of the graph, we introduce the matrices of the 
"right" and "left" turns 

The spaces Co and Ci are spaces of functions that take values in the fundamental two- 
dimensional representation of the group PS'L(2,R). 

The "operators of the right and left turns," Rz and Lz, are 

Rz^RXz = y (6.3) 

Lz^LXz = (6.4) 



The operator T{u) acts on the space Ci as follows. Let cr (respectively, ei) be the 
oriented edge to which the oriented edge is naturally mapped by the right (respectively, 
left) turn when going along consecutive edges of a path. For m G C and e G L, we obtain 

T{u)vg ■ Cz = vg^ ■ Cr + vg^ ■ Cl, 

where 

vg^ = uRzVg and vg^^ = uLzVg. 

A geodesic is a closed primitive path in the graph T. To each such path we set into 
the correspondence the product of matrices Pzi- z„ = Lz„Lz„^j^Rz„^2 ' ' ' -^^2-^21 5 where the 
matrices Lz- or Rz- are inserted depending on which turn — left or right — the path is going 
on the corresponding step. 

The matrices L and R are torsion potentials from Definition |^. However, we can modify 
the original graph T at each vertex disconnecting edges at the vertex and connecting them 

^ Note that 7ri(5) is isomorphic to 7ri(r). 
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by three "short" edges forming a triangle subsequently erasing one (any) of the new edges 
in order to preserve the number of loops of the graph (see Fig. 6). 




Fig. 6. "Blowing up" vertices of the initial graph. 



The matrices L and R then become the potential matrices corresponding to the "short" 
oriented edges of the resulting graph T. 

Proposition 4 [|14| There is a one-to-one correspondence between the set of all primitive 
closed paths {Pz^...z„} in the graph T and closed geodesies {7} on the Riemann surface. 
Moreover, the length ^(7) of a geodesic is determined by the relation 

G(7) = 2cosh(/(7)/2) = trP,,...,„. (6.5) 



Remark 4 Because both matrices (|6.3| ) and ( |6.4|) have the form | ^ ^ |, then any 

product of such matrices will be of the same form; meanwhile, in a diagonal term of any 
such product, all summands enter with the plus sign and, in an antidiagonal term, with 
the minus sign. Then, for closed geodesies around holes (round-the-face geodesies), we 
obtain /(7) = \ where the sum ranges all boundary edges of the face (with the 

proper multiplicities). 



Definition 12 The Ihara-Selberg L-function L{u) for the fat graph T is 

L(M) = ndet-^(/-M"P,,..,J (6.6) 

{^} 

(cf. (|3.1| )) where the product runs over all primitive closed paths, n being the length of a 
path measured in terms of the distance on the universal covering tree. 

Then, from Theorem |I|, we obtain that the Ihara-Selberg L-function ( |6.6| ) is a rational 
function in variables e ^'^^ and u, and 

L{u) = det -\I - T{u)) = det -^A{u). (6.7) 

The product ( |6.6| ) is absolutely convergent in the circle \u\ < 1/2 min{| e ^^'|}. How- 
ever, following Theorem |l], it possesses a unique analytic continuation into the whole C 
except a finite number of singular points. 
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6.1 Selberg trace formula and distribution of geodesies 

In the case of closed Riemann surface of constant negative curvature, strong results con- 
cerning the distribution of closed geodesies over lengths have been obtained (see ^7\, |IT]). 

Definition 13 The zeta-function for closed geodesies on the (punctured) surface S is 

C(.)=n(l-e-^'^'")-\ (6.8) 

n 

where /i is a constant that are related to the asymptotic distribution of geodesic lengths. 

If vr(r) is the number of closed (primitive) geodesies with the length at most T, then 
the following assertion holds true (see, e.g., |2^, [11[). 

Claim 2 There exists such constant h that 

lim — , „ ^ = 1. 

T^oo e'^^/hT 

The proof of Claim ^ follows from the ^-function analyticity properties. 

Proposition 5 [^] ({s) has an extension to C as a meromorphic function such that: 

(a) ({s) has no zeros or poles in 3?(s) > 1, s 7^ 1; 

(b) C(s) has a simple pole at s = 1. 

There are two alternative proofs of Proposition ^: one is based on the celebrated 
Selberg trace formula, which, in its simplest form, can be formulated as follows. 



Proposition 6 (for proofs, see p3|, |^). Given the Laplace-Beltrami operator A, which 



acts on the hyperbolic constant negative curvature metric space and is the linear second- 
order [unbounded) partial differential operator with the discrete spectrum = Aq < Ai < 
A2 < . . ., and the set of closed geodesies with lengths {/„} on the Riemann surface, the 
Selberg trace formula relates these two sequences of real numbers {A„} and {/„} by the 
formula 

Y.f{^\n-l/A) + / /dt = ^c.(/(Z„)) + 2^^^^ r rtanh(r)/(r)rfr, 

n n J -co 

where / : R — R «s a C°° function of compact support, f is its Fourier transform, and 

( f(] \\ — ^nfjkln) 

"^^^"^^"tiSinh(A;/„/2)' 
while Area {V)/27t = 2{g — 1) for a compact surface. 
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Alternative proof presented in [E^] just uses the methods of symbohc dynamics and can 



be rather close to the combinatorial graph description. Here we formulate two problems. 



Problem 3 To relate the graph description of zeta- (or L-)function (|3.1| ) of the Te- 
ichmiiller space with the standard zeta-function 



Problem 4 Few explicit formulas containing the geodesies in the open Riemann surface 
case are known. Worth mentioning is the paper where the formulas concerning sets 
of simple {i.e., non-self-intersecting) geodesies were obtained for a punctured torus and 
for a pair of pants. Are there generalizations of these formulas to a case of surfaces with 
holes? 



6.2 Classical projective (modular) transformations 

In [|^, the projective transformations on graphs that are the mapping class group (mod- 
ular) transformations were obtained. They correspond to natural operations called the 
flips, or Whitehead moves, which are elementary transitions between graphs (actually, 
between neighbor cells of the simplicial complex whose higher dimensional cells label 
combinatorial types of three- valent fat graphs of the given Riemannian genus g and n). 



The corresponding transformation of the variables is nonlinear (see M, IH] for the 



explicit expressions), but the geodesic length is a modular invariant function. 



Lemma 3 Atu = 1, the function L{u) is modular invariant, i.e., does not depend on the 
particular form of the representing graph. 



Note, however, that the graph length n of a geodesic varies under the modular trans- 
formations, so we are still unable to define a complete analogue of the zeta function (|6.8|) 
in terms of graphs. However, constructing L-functions ( |6.6D makes sense in the sym- 
bolic dynamics setting where one customarily introduce additional modular-noninvariant 
partition of a Riemann surface, which would corresponds to a graph decomposition (see 
Chaps. 5 and 6 in |TT|). 

Remark 5 As for the function L~^{u), the determinant in ( |6.6|) is a Laurent polynomial 
of no more than second order in e^^'^^ for each Zi. Moreover, from Lemma ^ it follows 
that the modular-invariant expression for L(l) can depend only on the modular invariants 
"perimeters" of holes. This imposes severe restrictions on a possible form of L^^(l). Say, 
for the moduli spaces i of the Riemann surfaces with one puncture (hole), there exists 
only one such parameter, pi. Taking into account the global symmetry Xi — Xj, we find 
that 

Lg^{l) = acosh(pi) + 6cosh(pi/2) + c, a + b + c = 0, 
where only the coefficients a, b, and c depend on the genus of the surface. 
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Example 4 1. For the moduli space Mi^i (the torus with one puncture), there are three 
parameters x, y, and z and one modular invariant pi = x + y + z, 

= -8(cosh(x + y + z)-l). 

2. The moduli space Mo,3 (the sphere with three punctures). There are three param- 
eters X, y, and z, and three restrictions x + y = ps, x + z = p2, y + z = pi. Then 

= -2= (co* (i±l) - l) (co* (i±i) - l) (cosh (i±i) - 1) . 
7 Conclusions 

Our main result is the relation ([4.2CI|) connecting the L-functions with the scattering data. 
It provides an analogue of the Selberg trace formula for the discrete non-compact case of 
the graphs. 

Spectral properties of the S'-matrix studied in Sec. 4 deserve further investigation, 
especially for the case of unitary and Schrodinger type potentials. 

One can hope to apply the technique of this paper to the continuous non-compact case 
in order to obtain analogous formulas for the continuum scattering processes. 

Eventually, some r-matrix structure must be hidden in these scattering processes. It 
is interesting how it can be interpreted from the standpoint of the integrable models. 

After this paper has been completed we were aware of paper [|^, where the similar 
results concerning magnetic fluxes on graphs have been obtained within the functional 
integral standpoint. 
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